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Abstract 

A search for neutral Higgs bosons of the IVIinimal Supersymmetric Standard IVIodel (IVISSIVI) is 
reported. The analysis is based on a sample of proton-proton collisions at a centre-of-mass energy 
of 7 TeV recorded with the ATLAS detector at the Large Hadron Collider. The data were recorded 
in 201 1 and correspond to an integrated luminosity of 4.7 fb~^ to 4.8 fb~^ Higgs boson decays into 
oppositely-charged muon or r lepton pairs are considered for final states requiring either the presence 
or absence of 6-jets. No statistically significant excess over the expected background is observed and 
exclusion limits at the 95% confidence level are derived. The exclusion limits are for the production 
cross-section of a generic neutral Higgs boson, (p, as a function of the Higgs boson mass and for 
h/A/H production in the MSSM as a function of the parameters and tan^ in the m™'''' scenario 
for ruA in the range of 90 GeV to 500 GeV. 
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Abstract: A search for neutral Higgs bosons of the Minimal Supersymmetric Standard 
Model (MSSM) is reported. The analysis is based on a sample of proton-proton collisions 
at a centre-of-mass energy of 7 TeV recorded with the ATLAS detector at the Large Hadron 
Collider. The data were recorded in 2011 and correspond to an integrated luminosity of 
4.7 fb~^ to 4.8 fb~^. Higgs boson decays into oppositely-charged muon or r lepton pairs are 
considered for final states requiring either the presence or absence of 6-jets. No statistically 
significant excess over the expected background is observed and exclusion limits at the 95% 
confidence level are derived. The exclusion limits are for the production cross-section of 
a generic neutral Higgs boson, (p, as a function of the Higgs boson mass and for h/A/H 
production in the MSSM as a function of the parameters and tan f3 in the m™*^^ scenario 
for ruA in the range of 90 GeV to 500 GeV. 
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1 Introduction 

Discovering the mechanism responsible for electroweak symmetry breaking is one of the 
major goals of the physics programme at the Large Hadron Collider (LHC) [1]. In the 
Standard Model this mechanism requires the existence of a single scalar particle, the Higgs 
boson [2-6]. The recent discovery of a particle compatible with the Higgs boson at the 
LHC [7, 8] provides further evidence in support of this simple picture. Even if this recently 
discovered particle is shown to have properties very close to the Standard Model Higgs 
boson, there are still a number of problems that are not addressed. For instance, quantum 
corrections to the mass of the Higgs boson contain quadratic divergences. This problem 
can be solved by introducing supersymmetry, a symmetry between fermions and bosons, 
by which the divergent corrections to the Higgs boson mass are cancelled. 

In the Minimal Super symmetric Standard Model (MSSM) [9, 10], two Higgs doublets 
are necessary, coupling separately to up-type and down-type fermions. This results in five 
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physical Higgs bosons, two of which are neutral and CP-even [h, -ff)^, one of which is 
neutral and CP-odd (A) , and two of which are charged {H^ ) . At tree level their properties 
can be described in terms of two parameters, typically chosen to be the mass of the CP- 
odd Higgs boson, tha, and the ratio of the vacuum expectation values of the two Higgs 
doublets, tan/3. In the MSSM, the Higgs boson couplings to r leptons and 6-quarks are 
strongly enhanced for a large part of the parameter space. This is especially true for large 
values of tan /3, in which case the decay of a Higgs boson to a pair of r leptons or 6-quarks 
and its production in association with 6-quarks play a much more important role than in 
the Standard Model. 

The results presented in this paper are interpreted in the context of the m™^^ bench- 
mark scenario [11]. In the scenario the parameters of the model are chosen such 
that the mass of the lightest CP-even Higgs boson, h, is maximised for a given point in 
the m^~tan/3 plane, under certain assumptions. This guarantees conservative exclusion 
bounds from the LEP experiments. The sign of the Higgs sector bilinear coupling, fi, is 
generally not constrained, but for the analysis presented in this paper /i > is chosen as 
this is favoured by the measurements of the anomalous magnetic dipole moment of the 
muon [12]. 

The masses of the Higgs bosons in this scenario are such that for large values of tan /3 
two of the three neutral Higgs bosons are closely degenerate in mass: for niA ^ 130 GeV, 
nih — mA and mn ~ 130 GeV, whereas for iua ^ 130 GeV, mn — mriA and ruh — 130 GeV. 

The most common MSSM neutral Higgs boson production mechanisms at a hadron 
collider are the 6-quark associated production and gluon-fusion processes, the latter pro- 
cess proceeding primarily through a &-quark loop for intermediate and high tan /?. Both 
processes have cross-sections that increase with tan /?, with the 6- associated production 
process becoming dominant at high tan /3 values. The most common decay modes at high 
tan/3 are to a pair of 5-quarks or r leptons, with branching ratios close to 90% and 10%, 
respectively, across the mass range considered. The direct decay into two muons occurs 
rarely, with a branching ratio around 0.04%, but offers a clean signature. 

Previous searches for neutral MSSM Higgs bosons have been performed at LEP [13], 
the Tevatron [14-16] and the LHC [7, 8]. The recently observed Higgs-boson-like particle 
at the LHC [17, 18] is consistent with both the Standard Model and the lightest CP-even 
MSSM Higgs boson [19, 20]. In this paper a search for neutral MSSM Higgs bosons using 
4.7 fb~^ to 4.8 fb~^ of proton-proton collision data collected with the ATLAS detector [21] 
in 2011 at the centre-of-mass energy of 7TeV is presented. The and t~^t~ decay 

modes are considered, with the latter divided into separate search channels according to 
the subsequent r lepton decay modes. Events from each channel are further classified 
according to the presence or the absence of an identified 6-jet. 



^By convention the lighter CP-even Higgs boson is called h, the heavier CP-even Higgs boson is called 

H. 
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2 The ATLAS detector 



The ATLAS experiment at the LHC is a multi-purpose particle detector with a forward- 
backward symmetric cylindrical geometry and nearly 47r coverage in solid angle [21]. It 
consists of an inner tracking detector surrounded by a thin superconducting solenoid pro- 
viding a 2 T axial magnetic field, electromagnetic and hadronic calorimeters, and a muon 
spectrometer. The inner tracking detector covers the pseudorapidity range \r]\ < 2.5^. 
It consists of silicon pixel, semi-conductor micro-strip, and transition radiation tracking 
detectors. Lead/liquid- argon (LAr) sampling calorimeters provide electromagnetic (EM) 
energy measurements with high granularity. A hadronic (iron/scintillator-tile) calorimeter 
covers the central pseudorapidity range (|?7| < 1.7). The end-cap and forward regions are 
instrumented with LAr calorimeters for both EM and hadronic energy measurements up to 
|?7| = 4.9. The muon spectrometer surrounds the calorimeters and incorporates three large 
air-core toroid superconducting magnets with bending power between 2.0 Tm and 7.5 Tm, 
a system of precision tracking chambers and fast detectors for triggering. A three-level 
trigger system is used to select events. The first-level trigger is implemented in hardware 
and uses a subset of the detector information to reduce the rate to at most 75 kHz. This is 
followed by two software-based trigger levels that together reduce the event rate to approx- 
imately 300 Hz. The trigger requirements were adjusted to changing data-taking conditions 
during 2011. 

3 Data and Monte Carlo simulation samples 

The data used in this search were recorded by the ATLAS experiment during the 2011 LHC 
run with proton-proton collisions at a centre-of-mass energy of 7TeV. They correspond 
to an integrated luminosity of 4.7 fb~^ (r"*"T~ channels) or 4.8 fb~^ {fi'^fi~ channel) after 
imposing the data quality selection criteria to require that all relevant detector sub-systems 
used in these analyses were operational. 

Higgs boson production: The Higgs boson production mechanisms considered are 
gluon-fusion and 6-associated production. The cross-sections for the first process have been 
calculated using HIGLU [22] and ggh@nnlo [23]. For 6-associated production, a matching 
scheme described in reference [24] is used to combine four- flavour [25, 26] calculations and 
the five-flavour bbh@nlo [27] calculation. The masses, couplings and branching ratios of 
the Higgs bosons are computed with FeynHiggs [28]. Details of the calculations and the 
associated uncertainties due to the choice of the value of the strong coupling constant, the 
parton distribution function and the factorisation and renormalisation scales can be found 
in reference [29]. Gluon-fusion production is simulated with POWHEG [30], while 6-quark 
associated production is simulated with SHERPA [31]. 

^ ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in 
the centre of the detector and the z-axis along the beam pipe. The x-axis points from the IP to the centre 
of the LHC ring, and the j/-axis points upwards. Cylindrical coordinates (r, <j}) are used in the transverse 
plane, (p being the azimuthal angle around the beam pipe. The pseudorapidity is defined in terms of the 
polar angle as rj = ~ lntan(S/2). 
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The h/A/H — t- /-i^/i^ and h/A/H — )• t^t~ modes are considered for the decay of the 
Higgs boson. The A boson samples are generated for both production mechanisms and 
are also employed for modelling H and h production. The differences between CP-even 
and CP-odd eigenstates are negligible for this analysis. The signal modelling for a given 
combination of itla and tan (3 takes into account all three Higgs bosons /i, H and A, by 
adding the A boson mass samples corresponding to m/j, uih and ttia according to their 
production cross-sections and masses in the m™^^, /i > MSSM benchmark scenario. 

For the t^t~ decay mode, 15 Monte Carlo samples are generated with Higgs boson 
masses in the range of 90GeV to 500 GeV and tan^ = 20. These are scaled to the appro- 
priate cross-sections for other tan j3 values. The simulated signal samples with ttia closest 
to the computed mass of H and h are used for the H and h bosons. The increase in the 
Higgs boson natural width with tan /3, of the order of 1 GeV in the range considered, is 
negligible compared to the experimental mass resolution in this channel, which is always 
above 10 GeV. 

For the decay mode, seven samples are generated with Higgs boson masses in the 

range of 110 GeV to 300 GeV and tan/3 = 40. Additionally, to study the tan/3 dependence 
of the width of the resonance, signal samples are generated for both production modes for 
rriA = 150 GeV and 250 GeV, each at tan /3 = 20 and tan /3 = 60. Since the mass resolution 
is better in this channel, signal distributions are obtained using an interpolation procedure 
for different intermediate m^-tan/3 values, as described in section 5. 

The generated Monte Carlo samples for the h/A/H — t- r'^r" decay modes are passed 
through the full GEANT4 [32, 33] detector simulation, while the samples for the h/A/H — )• 
decay mode are passed through the full GEANT4 detector simulation or the "fast" 
simulation, ATLFAST-H [32], of the ATLAS detector. 

Background processes: The production of W and Z/^* bosons in association with jets 
is simulated using the ALPGEN [34] and PYTHIA [35] generators. PYTHIA is also used 
for the simulation of hh production, but through an interface, which ensures that the simu- 
lation is in agreement with 6-quark production data [36, 37]. The if production process is 
generated with MC@NLO [38] and POWHEG [39, 40]. MC@NLO is used for the genera- 
tion of electroweak diboson (W^VF, WZ^ ZZ) samples. Single-top production through the 
s- and t-channels, and in association with W bosons, is generated using AcerMC [41]. For 
all event samples described above, parton showers and hadronisation are simulated with 
HERWIG [42] and the activity of the underlying event with JIMMY [43]. The loop-induced 
gg —7- WW processes are generated using gg2WW [44]. The following parton distribution 
function sets are used: CTIO [45] for MC@NLO, CTEQ6L1 [46] for ALPGEN and modified 
leading-order MRST2007 [47] for PYTHIA samples. 

Decays of r leptons are simulated using either SHERPA or TAUOLA [48]. Initial- 
state and final-state radiation of photons is simulated using either PHOTOS [49] or, for 
the samples generated with SHERPA, PHOTONS-F-F, which is a part of SHERPA. The 
Z/7* — )• r+T~ background processes are modelled with a r-embedded Z/^* — )• data 
sample described in section 6. All generated Monte Carlo background samples are passed 
through the full GEANT4 simulation of the ATLAS detector. 
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The signal and background samples are reconstructed with the same software as used 
for data. To take into account the presence of multiple interactions occurring in the same 
and neighbouring bunch crossings (referred to as pile- up), simulated minimum bias events 
are added to the hard process in each generated event. Prior to the analysis, simulated 
events are re-weighted in order to match the distribution of the number of pile-up interac- 
tions per bunch crossing in the data. 

4 Physics object reconstruction 

An electron candidate is formed from energy deposits in the electromagnetic calorimeter 
associated with a charged particle track measured in the inner detector. Electrons are 
selected if they have a transverse energy > 15GeV, lie within \r]\ < 2.47, but outside 
the transition region between the barrel and end-cap calorimeters (1.37 < \r]\ < 1.52), and 
meet quality requirements based on the expected shower shape [50] . 

A muon candidate is formed from a high-quality track measured in the inner detector 
matched to hits in the muon spectrometer [51]. Muons are required to have pT > lOGeV 
and to lie within |?/| < 2.5. In addition, the point of closest approach of the inner detector 
track must be no further than 1 cm from the primary vertex^, as measured along the z-axis. 
This requirement reduces the contamination due to cosmic-ray muons and beam-induced 
backgrounds. 

Identified electrons and muons are isolated if there is little additional activity in the 
inner detector and the calorimeter around the electron or muon. The scalar sum of the 
transverse momenta of all tracks from the same vertex as the lepton, with above 1 GeV 
and located within a cone with radius parameter^ AR = 0.4 around the lepton direction, 
must be less than 6% of the lepton momentum for the r+r" channels, or less than 10% 
for the /^"^/U" channels. The sum excludes the track associated to the lepton itself. In 
addition, a pile-up-dependent requirement on the calorimetric energy within a cone of 
radius AR = 0.2 around the lepton direction, excluding the energy associated to the 
lepton itself, is applied in the t~^t~ channels using a calculation of the median transverse 
energy density from Iow-Et jets [52, 53]. 

Reconstructed electrons and muons are used to identify the leptonic decays of r leptons 
(riep). The decay of a r lepton into an electron (muon) and neutrinos is denoted by Te (r^). 

Jets are reconstructed using the anti-A;f algorithm [54] with a radius parameter R = 0.4, 
taking three-dimensional noise-suppressed clusters of calorimeter-cell energy deposits [55] 
as input. The jet energy is calibrated using a combination of test-beam results, simulation 
and in situ measurements [56]. Jets must satisfy Ej- > 20 GeV and \rj\ < 4.5. Rare events 
containing a jet with associated out-of-time activity or calorimeter noise are discarded. To 
reduce the effect of pile-up, at least 75% of a jet's transverse momentum, as measured by 
the scalar sum of the transverse momenta of the associated tracks, must come from tracks 
matched to the primary vertex. A multivariate algorithm based on a neural network is 

^The primary vertex is defined as the vertex with the largest Epx of the associated tracks. 

*Ai? — \J (A77)^ + (A(/))^, where Ary is the difference in pseudorapidity of the two objects in question 
and A0 is the difference between their azimuthal angles. 
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used in this analysis to tag jets, reconstructed within \r]\ < 2.5, that are associated with 
the hadronisation of 6-quarks. The neural network makes use of the impact parameter of 
associated tracks and the reconstruction of b- and c-hadron decay vertices inside the jet [57]. 
The 6-jet identification has an efficiency of about 70% in tt events, unless otherwise stated. 
The corresponding rejection factors are about 5 for jets containing charm hadrons and 
about 130 for light-quark or gluon jets. 

Hadronic decays of r leptons (rhad) are characterised by the presence of one, three, or 
in rare cases, five charged hadrons accompanied by a neutrino and possibly neutral hadrons, 
resulting in a collimated shower profile in the calorimeters with only a few nearby tracks. 
The visible r decay products are reconstructed in the same way as jets, but are calibrated 
separately to account for the different calorimeter response compared to jets. Information 
on the collimation, isolation, and shower profile is combined into a boosted-decision-tree 
discriminant to reject backgrounds from jets [58]. In this analysis, three selections are 
used — "loose", "medium" and "tight" — with identification efficiency of about 60%, 45% 
and 35%, respectively. The rejection factor against jets varies from about 20 for the loose 
selection to about 300 for the tight selection. A Thad candidate must lie within |r/| < 2.5, 
have a transverse momentum greater than 20 GeV, one or three associated tracks (with 
Pt > IGeV), and a total charge of ±1. Dedicated electron and muon veto algorithms are 
used. 

When different objects selected according to the above criteria overlap with each other 
geometrically (within AR = 0.2), only one of them is considered for further analysis. The 
overlap is resolved by selecting muon candidates, electron candidates, Thad candidates and 
jet candidates in this order of priority. 

The magnitude and direction of the missing transverse momentum, E^^^^, is recon- 
structed including contributions from muons and energy deposits in the calorimeters [59]. 
Clusters of calorimeter-cell energy deposits belonging to jets, Thad candidates, electrons, 
and photons, as well as cells that are not associated to any object, are treated separately 
in the E^^^^ calculation. The contributions of muons to EJ^^^^ are calculated differently 
for isolated and non-isolated muons, to account for the energy deposited by muons in the 
calorimeters. 

5 The fJ.~^n~ decay channel 

Signal topology and event selection: The signature of the h/A/H — t- ^'^ ^~ decay 
is a pair of isolated muons with high transverse momenta and opposite charge. In the 
6-associated production mode, the final state can be further characterised by the presence 
of one or two Iow-Et 6-jets. The missing transverse momentum is expected to be small and 
of the order of the resolution of the EJ^^^^ measurement. The n~^lJ-~ decay channel search 
is complicated by a small branching ratio and considerable background rates. 

Events considered in the analysis must pass a single-muon trigger with a trans- 

verse momentum threshold of 18 GeV. At least one reconstructed muon is required to be 
matched to the T]-(p region of the trigger object and to have pT > 20 GeV. At least one 
additional muon of opposite charge and with > 15 GeV is required. A muon pair is 
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Figure 1. The invariant mass distribution of the two muons of the h/A/H — >■ jX^ ^~ search for 
the 6-tagged (left-hand side) and the 6- vetoed samples (right-hand side). The data are compared 
to the background expectation and a hypothetical MSSM signal (m^ = 150 GeV, tan/3 = 40). 
Simulated backgrounds are shown for illustration purposes. The background uncertainties shown 
here are statistical in nature due to the finite number of simulated events. The contributions of 
the backgrounds Z/^* — > e+e^, W + jets and those of all diboson production processes but WW 
production are combined and labelled "Other electroweak" . 



formed using the two highest-px muons of opposite charge. This muon pair is required 
to have an invariant mass greater than 70 GeV. In addition, events are required to have 
^miss ^ QgY niuons considered here must be isolated, as defined in section 4. 

The large background due to Z/^* production can be reduced by requiring that the 
event contains at least one identified 6-jet. Events satisfying this requirement are included 
in the 6-tagged sample, whereas events failing it are included in the 6-vetoed sample. The 
/i"*"/!" invariant mass distribution, ?n^^, is shown separately for the 6-tagged and the b- 
vetoed samples in figure 1. For illustration purposes only, the distributions of simulated 
backgrounds and an assumed MSSM neutral Higgs boson signal with ruA = 150 GeV and 
tan^ = 40 are shown in the same figure. A hypothetical signal would be present as 
narrow peaks on top of the high-mass tail of the Z boson superimposed on a continuous 
contribution from non-resonant backgrounds such as ti. The Z/^* process contributes 
to the total background with a relative fraction of about 99% (51%) for the 6-vetoed (6- 
tagged) sample for events in the m^^ range of 110 GeV to 300 GeV, which is most relevant 
to the Higgs boson searches in this channel. In the 6-vetoed sample the remaining non- 
resonant background is composed of tt, W~^W~ and hh events while M events dominate the 
non-resonant background in the 6-tagged sample. 

Background modelling: The background in the channel is estimated from data. 

By scanning over the /U"*"/!" invariant mass distribution, local sideband fits provide the ex- 
pected background estimate in the mass region of interest. To this end, a parameterisation 
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of the background shape is fitted to the /^^/U^ invariant mass distribution. Search windows 
are defined around each of the neutral Higgs bosons and are excluded from the fit. This 
results in one or two windows due to the mass degeneracy among the Higgs bosons. The 
widths of the search windows are motivated by the expected signal width for each point 
in the scanned m^i-tan (3 grid and account for asymmetries in the signal invariant mass 
distribution. The upper and lower boundaries of the search windows are defined by the 
m^^ values where the cross-section predictions of the signal model are 10% of their max- 
imum. The lower and upper outer boundaries of the sidebands vary between 98-118 GeV 
and 160-400 GeV, respectively, depending on rriA and m/,. 

The parameterisation of the invariant mass distribution, /^(x), is given by 



Jb {x \Nb,A, B, mz,rz, a) = Nb ■ [fz {x \A, B, mzJz) ® -^G [x |0, g)\ , 



(5.1) 



where x is the invariant mass, ® the convolution operator and J-g (x |0, cr) the Gaussian 
distribution with variable x, mean and variance . The function describing the Zj^* 
production is 



1 x^ — mi x^ 
~^ + To" 



h{x\A,B,mz,Tz) = A^^B- ^ + - ^2 • (5-2) 



This is convolved with the Gaussian distribution accounting for the finite mass resolution. 
The function fz is a simplification of the pure 7* and Z propagators, including Z-j* 
interference contributing to the process qq —7-^/7* — t- and hence in principle only 

describes the background from Z/^* production. The parameterisation /b is found to be 
a good approximation of the shape of the total /i"*"/^" background even in the 6-tagged 
sample, which has non- negligible contributions from physics processes other than Z/'j* 
production. 

In total, the fit function, /b, has six parameters. The natural width of the Z boson, Tz, 
is fixed to Tz = 2.50 GeV, whereas the remaining parameters are unconstrained. Parameter 
Nb describes the total normalisation of the curve and parameters A and B represent the 
relative normalisations of the 7* and Z--y* contributions with respect to the Z term. 
Finally, mz represents the mass of the Z boson and the parameter a specifies the mean 
pair mass resolution. 

For every point on the m^-tan/3 grid, a binned likelihood fit of Jb to the data is 
performed to estimate the five unconstrained parameters and consequently the total back- 
ground estimate. 

The background model is validated with x^-based goodness-of-fit studies. In addition, 
the background model is extended by polynomials of different orders to test if additional 
degrees of freedom change the goodness of the fit, which would hint at problems in the 
shape modelling. Further validation of the capability of the model to describe the shape 
of the data is performed by varying the fit ranges for certain mass points and accounting 
for the fit residuals. The goodness-of-fit studies confirm a good background modelling for 
both the 6-vetoed and the 6-tagged sample. The uncertainty on the background estimate 
is obtained from a variation of the fitted background function within its 68% confidence 



8 



level (CL) uncertainty band. This results in an uncertainty of 5% (2%) on the expected 
background yield for the 6-tagged (6-vetoed) sample. 

Signal modelling: The h/A/H — t- /i signal is expected to appear as narrow peaks 
in the invariant mass distribution, as depicted in figure 1. The resolution in the 

relevant mass range is typically 2.5% to 3%, and numerous mass points are needed for 
a complete mass scan. In addition, the influence of tan/3 on the reconstructed width of 
the signal invariant mass distribution needs to be taken into account. The natural widths 
of the MSSM neutral Higgs bosons increase with tan /3. The reconstructed width can be 
sensitive to this variation because of the good experimental mass resolution. 

To interpolate between the different signal samples obtained from a limited number of 
simulated signal masses, the signal ^i'^ ^~ invariant mass distribution is parameterised with 



where x represents the invariant mass. The parameterisation consists of a Breit- 

Wigner function describing the signal peak convolved with a Gaussian distribution, J^q, 
accounting for the finite mass resolution and a Landau function, with low-mass tail 
which models the asymmetric part of the signal invariant mass distribution. 

The function fs is characterised by six parameters. The width of the Breit-Wigner 
function, F, is fixed to the theoretical predictions calculated with FeynHiggs [28]. The 
remaining five parameters are unconstrained. The overall normalisation parameter is Ns 
and c specifies the relative normalisation of the Landau function with respect to the Breit- 
Wigner function. Parameter m specifies the mean of the Breit-Wigner and the Landau 
distributions, a determines the width of the Gaussian distribution and ? represents the 
scale parameter of the Landau function. 

The function fs is fitted to each signal sample available from simulation. Each fit 
results in a set of fitted parameters, (iV5, m, tr, c, depending on the point in the niA- 
tan /3 plane. The dependence of this set on and tan /3 is parameterised with polynomials 
of different orders. The resulting polynomials provide a set of parameters which in addition 
to the predicted natural width, F, fully define the normalised probability density function 
for an arbitrary point in the m^-tan /3 plane. 

This procedure is used to generate invariant mass distributions for signal masses from 
120 GeV to 150 GeV in 5 GeV steps and from 150 GeV to 300 GeV in 10 GeV steps, as weh 
as for tan/3 values from 5 to 70 in steps of 3 or 5. Higgs boson masses below 120 GeV were 
not considered because the background model does not provide precise estimates in the 
mass region close to the Z boson peak. For both the 6-tagged and the 6-vetoed samples 
the interpolated and normalised probability density functions are obtained separately for 
the Higgs boson production from gluon-fusion and in association with 6-quarks. As for 
the background, the uncertainty on the signal prediction from the fit is obtained from its 
68% CL uncertainty band. The resulting uncertainty is estimated to be 10% to 20% of the 
signal event yield. 



fs {x \Ns, m, F, a, c, = Ns 



1 



(g) Jb {x |0,cr) 




+ c • Jl (-X |m,?) , 



(5.3) 
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Figure 2. Invariant mass distribution of data and predicted background from sideband fits to the 
data shown for the signal mass point at ttia — 150 GeV and tan/3 = 40 for the 6-tagged (left-hand 
side) and the 6-vetoed samples (right-hand side) of the h/A/H — > iJ^fi^ final state. The ratio of the 
data to the predicted background, labelled D/B, and the bin- by-bin significances of the deviations 
of the data from the background prediction, labelled a, are shown beneath. 



Results: Figure 2 compares the data with the background estimate predicted from side- 
band fits in both the 6-tagged and 6- vetoed samples for the signal mass point rriA = 150 GeV 
and tan /3 = 40. The data fluctuate around the background prediction leading to local bin- 
by-bin significances that are typically less than 2 a. Table 1 shows the number of observed 
events in the fit range around the mass point rriA = 150 GeV compared to the number 
of background events predicted by the sideband fits. The observed numbers of events are 
compatible with the expected yield from Standard Model processes within the uncertain- 
ties. 



6 The T+T decay channel 

The h/A/H — t- t'^t~ decay mode is analysed in several categories according to the r lepton 
decay final-state combinations. The four decay modes considered here are: TeT^ (6%), 
TeThad (23%), T^Thad (23%) and Thad''"had (42%), where percentages in the parentheses denote 
the corresponding branching ratios. The combination of rgThad and r^Thad is referred to as 

T'lcpThad- 

6.1 Common background estimation and mass reconstruction techniques 

T-embedded Z/'^* — >■ data: Z/'^* — )■ t'^t~ events form a largely irreducible 

background to the Higgs boson signal in all final states. It is not possible to select a 
Z/^* — )• T^T~ control sample which is Higgs boson signal-free. However, Z/7* — )• /U"*^^" 
events can be selected in data with high purity and without significant signal contamination. 
Furthermore, the event topology and kinematics are, apart from the r lepton decays and 
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6-tagged sample 


6-vetoed sample 


Mass Point 




rriA = 


= 150 GeV 


Fit Range 




110- 


-200 GeV 


Background 




980 ± 50 


35900 ± 600 


Signal rriA = 


ioU LjeV , tan p = 4U 






bb{h/A/H^ 




28 ± 2 tl 


271 ± 22 til 


gg ^h/A/H- 




2.3 ±0.3 ±0.4 


141 ± 10 til 


Data 




985 


36044 



Table 1. Observed number of data events and the expected number of signal and background 
events in the h/A/H — > ^i'^ channel for one of the considered signal mass points. The number of 
background events is predicted from sideband fits to the /i"*"/!" invariant mass distribution in the 
fit range around the signal mass point niA ~ 150 GeV for both the h/A/H — >• 6-vetoed and 

6-tagged samples. The number of expected signal events produced in gluon-fusion or in association 
with 6-quarks is shown for tan/3 = 40. The quoted uncertainty for the background estimate is the 
statistical uncertainty obtained from the fit. For the signal estimate, the uncertainty from the fit 
is quoted first and then separately the uncertainty from other sources discussed in section 7. 

the different masses of r leptons and muons, identical to those of Z/"f* — >■ t'^t" events. 
Therefore Z/^* — )■ /-i^/i^ events are selected in data and modified using a r-embedding 
technique, in which muons are replaced by simulated r leptons. The hits of the muon tracks 
and the associated calorimeter cells in a cone with radius parameter Ai? = 0.1 around the 
muon direction are removed from the data event and replaced by the detector response from 
a simulated Z/^* — )■ t~^t~ event with the same kinematics. The event reconstruction is 
performed on the resulting hybrid event. Only the r decays and their detector response are 
taken from the simulation, whereas the underlying event kinematics and the associated jets 
are taken from the data event. The procedure treats consistently the effect of r polarisation 
and spin correlations. The event yield of the embedded sample after the selection of the 
r decay products is normalised to the corresponding event yield obtained in a simulated 
Z/7* — >• r+T~ sample. This procedure has been validated as described in references [7, 60]. 

Systematic uncertainties on the normalisation and shape of the embedded sample are 
derived by propagating variations of the Z/^* — )■ /U"^^" event selection and the muon 
energy subtraction procedure through the r-embedding process. Additional uncertainties 
are assigned due to the use of the Z/^* — )■ r'^r" cross-section and Monte Carlo acceptance 
prediction in determining the r-embedded — )• /U sample normalisation. These 
theoretical uncertainties are described in section 7. 

Jets misidentified as hadronic r decays: A fraction of jets originating from quarks 
or gluons are misidentified as Thad candidates. It has been shown in reference [61] that this 
misidentification fraction is higher in simulated samples than in data. To account for this 
difference, the Monte Carlo background estimate is corrected based on control samples. 
Details are presented for each decay channel separately. 
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The ABCD background estimation method: The estimation of the background 
from multi-jet processes is done from data using the ABCD method for all t~^t~ channels. 
Two uncorrelated variables are chosen to define four data regions, named A, B, C and D, 
such that one variable separates A and B from C and D, while the other separates A and 
C from B and D. The signal region is labelled A, and the other regions are dominated by 
background from multi-jet processes. An estimate of the background from these processes 
in the signal region, ua-, is: 

nA = nB -X = nB X rem, (6-1) 

no 

where n^, nc and n£, denote the populations of regions B, C and D, respectively. The 
populations of the B, C and D regions may need to be corrected by subtracting the 
estimated number of events that come from processes other than multi-jet production. 
This estimate is generally obtained from simulation. 

T+T" mass reconstruction: The invariant mass of the t~^t~ pair cannot be recon- 
structed directly due to the presence of neutrinos from the r lepton decays. Therefore, a 
technique known as the Missing Mass Calculator (MMC) is used to reconstruct the Higgs 
boson candidate mass [62] . This algorithm assumes that the missing transverse momentum 
is due entirely to the neutrinos, and performs a scan over the angles between the neutrinos 
and the visible r lepton decay products. For leptonic r decays, the scan also includes the 
invariant mass of the two neutrinos. At each point, the t~^t~ invariant mass is calculated, 
and the most likely value is chosen by weighting each solution according to probability 
density functions that are derived from simulated r lepton decays. This method provides 
a 13% to 20% resolution in the invariant mass, with an efficiency of 99% for the scan to 
find a solution. 

6.2 The h/A/H TeT^ decay channel 

Signal topology and event selection: Events in this channel must satisfy either a 
single-electron, single-muon or combined electron-muon trigger. The single-lepton triggers 
have pt thresholds of 20 GeV or 22 GeV for electrons, depending on the run period, and 
18 GeV for muons, while the combined trigger has a threshold of 10 GeV for the electron 
and 6 GeV for the muon. Exactly one isolated electron and one isolated muon of opposite 
electric charge are required, and the lepton pair must have an invariant mass exceeding 
30 GeV. The thresholds are 15 GeV for electrons and 10 GeV for muons in cases where 
the event is selected by the combined electron-muon trigger. These thresholds are raised 
to 24 GeV for electrons and 20 GeV for muons in cases where the event is selected by a 
single-lepton trigger only. 

The event sample is then split according to its jet flavour content. Events containing 
exactly one identified 6-jet are included in the 6-tagged sample. This is based on a 6-jet 
identification criterion with 75% efficiency in ti events. Events without an identified 6-jet 
are included in the 6-vetoed sample. The scalar sum of the lepton transverse momenta 
and missing transverse momentum is required to fulfil the following condition to reduce 
top quark and diboson backgrounds: E^^^^ + PT,e + Pt,/x < 125 GeV (< 150 GeV) for the 
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6-tagged (6-vetoed) sample. In order to further suppress these backgrounds and W — )• 
iv events, the opening angle between the two lepton candidates in the transverse plane 
must satisfy the condition Acpe^ > 2.0 (> 1.6) for the 6-tagged (6-vetoed) sample. In 
addition, the combination of the transverse opening angles between the lepton directions 
and the direction of E^^^^ is required to satisfy the condition J2e=e ^ A</)£;miss ^ > —0.2 
(> —0.4) for the 6-tagged (6-vetoed) sample. Finally, the scalar sum of the transverse 
energies of all jets, H^, is restricted to be below 100 GeV in the 6-tagged sample to further 
suppress backgrounds containing a higher multiplicity of jets, or jets with higher transverse 
momenta, than expected from the signal processes. Jets with |7y| < 4.5 and Et > 20 GeV 
are used to calculate the value of Ht. 

Estimation of the Z/^* — >■ t~^t~ background: The Z/j* — )■ t~^t~ background is es- 
timated by using the r-embedded Z/j* — )■ /U"*"^" event sample outlined in section 6.1. The 
use of multiple triggers with different pT thresholds has an effect on the lepton transverse 
momentum spectra, which is accurately reproduced by the trigger simulation. However, in 
the r-embedded Z/'j* — )• data there is no simulation of the trigger response for the 

decay products of the r leptons. This has an impact on the MMC mass distribution in the 
r-embedded Z/^* — t- fJ-'^f^" data, which is comparable to the statistical uncertainty in the 
6-vetoed sample, and negligible in the 5-tagged sample. For this reason the trigger selection 
is emulated for the 6-vetoed sample such that the trigger effect is adequately described. 
This emulation is based on the px-dependent trigger efficiencies obtained from data. 

Estimation of the tt background: The contribution of tt production is extrapolated 
from control regions which have purities of 90% (6-tagged sample) and 96% (6-vetoed 
sample). The selection criteria for these control regions are identical to the respective signal 
regions with two exceptions: at least two identified 6-jets are required, and the selection 
-ffx < 100 GeV is not applied. The multi-jet contributions to these control regions are 
estimated from data with the ABCD method; the other non-tt contributions are taken 
from simulation. The uncertainty on the normalisation obtained in this manner is 15% 
(30%) in the 6-tagged (6-vetoed) sample, primarily due to uncertainties on the 6-tagging 
efficiency and jet energy scale. 

Estimation of the multi-jet background: The multi-jet background is estimated us- 
ing the ABCD method, by splitting the event sample into four regions according to the 
charge product of the e/z pair and the isolation requirements on the electron and muon. 
These requirements are summarised in table 2. 

The systematic uncertainty of this method has been estimated by considering the sta- 
bility of the ratio VQjj;) in regions where the isolation requirements are varied, or where only 
the muon is required to fail the isolation. The resulting uncertainty on the normalisation 
is 14% (23%) in the 6-tagged (6-vetoed) sample. 

Smaller backgrounds from + jets, Z/^* — )• e''~e~, Z/^* — ^"""/x", diboson, and 
single-top processes are estimated from simulation. 

Results: The number of observed rgr^ events in data, along with predicted event yields 
from background processes, is shown in table 3. The observed event yield is compatible 
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Region Charge correlation Lepton isolation requirement 

A (Signal Region) Opposite sign isolated 

B Same sign isolated 

C Opposite sign anti-isolated 

D Same sign anti-isolated 



Table 2. Control regions for the estimation of the multi-jet background for the h/A/H — > TgT^ 
and h/A/H — >• riopThad samples: events are categorised according to the charge product of the two 
r leptons and the lepton isolation requirement. In the h/A/H — > riopThad channel isolation refers 
to the isolation of the electron or muon and in the h/A/H — > Tgr^ channel both the electron and 
muon are required to be isolated or anti-isolated, respectively. 
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Figure 3. MMC mass distributions for the h/A/H TeT^, final state. The MMC mass, m^^'^, 
is shown for the 5-taggcd (left-hand side) and 6- vetoed samples (right-hand side). The data are 
compared to the background expectation and an added hypothetical MSSM signal (niA = 150 GeV 
and tan/3 = 20). The background uncertainties include statistical and systematic uncertainties. 
The contributions of the backgrounds Z/^* — >• e+e", Z/j* — >• /i"*"//" and W -I- jets are combined 
and labelled "Other electroweak" ; in the case of the 6-tagged samples the contributions of diboson 
production processes are included as well. Background contributions from top quarks are included 
in "Other electroweak" for the 6-vetoed sample. 



with the expected event yield from Standard Model processes within the uncertainties. The 
MMC mass distributions for these events are shown in figure 3. 

6.3 The h/A/H "HepTTiad decay channel 

Signal topology and event selection: Events in the h/A/H — )• riepThad channel are 
selected using a single-lepton trigger with transverse momentum thresholds of 20 GeV or 
22 GeV for electrons, depending on the run period, and 18 GeV for muons. Each event must 
contain one isolated electron with > 25 GeV or one isolated muon with > 20 GeV. 
Events containing additional electrons or muons with transverse momenta greater than 
15 GeV or 10 GeV, respectively, are rejected in order to obtain an orthogonal selection to 
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6-tag 


ged sample 


6-vetoed sample 


Z/7* T+T- 


109 


±12 


11000 


±1000 


Ty + jets 


1.2 


H-i.i 
-0.9 


111 


±23 




1.1 


±0.8 


196 


+22 
-23 


tt 


56 


4-1 1 

-9 


150 


+60 

-50 


Single top 


16 


-4 


35 


±5 


Diboson 


3.9 


±0.7 


470 


±50 


Multi-jet 


15 


±11 


980 


±230 


iotai 


ZUi 


+20 
-19 


13000 


±1000 


Signal myi = 150 GeV, tan (3 = 20 










bb{h/A/H^ tt) 


18 


+4 
-5 


270 


+40 

-50 


gg -^h/A/H^ tt 


2.3 


±0.8 


143 


+23 
-21 


Data 


181 




12947 





Table 3. The number of events observed in data and the expected number of signal and background 
events of the h/A/H — > TeT^ channel. Simulated event yields are normalised to the integrated 
luminosity of the data sample, 4.7 fb^^. The predicted signal event yields correspond to a parameter 
choice of mA = 150 GeV and tan/3 = 20 and include both the 6-associated and the gluon-fusion 
production processes. 

those used in the Tgr^ and /U^/ti^" channels. One Thad with a charge of opposite sign to 
the selected electron or muon is required. The Thad identification criterion in use is the 
one with medium efficiency as introduced in section 4. The transverse mass of the lepton 
and the missing transverse momentum, mx, is required to be less than 30 GeV, to reduce 
contamination from W + jets and tt background processes. Here the transverse mass is 
defined as: 

mx = ^ 2p^^ Ef^^ {I - cos A(P) , (6.2) 

where p}^^ denotes the transverse momentum of the electron or muon and Acf) the angle 
between p^^^ and 

After imposing these selection criteria, the resulting event sample is split into two 
categories depending on whether or not the highest-£'T jet with |r/| < 2.5 is identified as a 
6-jet. Events are included in the 5-tagged sample if the highest-^Er jet is identified as a 5-jet 
and its Et is in the range of 20 GeV to 50 GeV. Events are included in the 6-vetoed sample if 
the highest-£'T jet fails the 5-jet identification criterion and the event has E!f^^^ > 20 GeV. 

Estimation of the W + jets background: Vr±jets events that pass the event selection 
criteria up to the mx requirement consist primarily of events in which the selected lepton 
originates from the W decay and a jet is misidentified as a Thad- To ensure a proper 
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estimation of the jet-to- Thad misidentification rate, the W + jets background normaUsation 
is corrected using control regions with high purity in Vl^ + jets events defined by requiring 
high transverse mass: 70 < rriT < 110 GeV. Separate control regions are used for the 
TeThad ^iiid T^T^ad Samples, as the kinematic selections are different. The correction factors 
derived from these control regions are = 0.587 it 0.009 for the electron channel and 
fw ~ 0-541 lb 0.008 for the muon channel, where the quoted uncertainties are statistical. 
The relative systematic uncertainty is estimated to be 5% by varying the mx boundary 
definition of the control region. The correction factors have been derived separately for the 
6-tagged and fe-vetoed samples; the numbers are in agreement between the two cases, but 
for the 6-tagged sample the statistical uncertainty is 17%. This statistical uncertainty is 
considered as an additional systematic uncertainty in the 6-tagged sample analysis. 

Estimation of the Z/7* T+T~/e+e~//i+/x~ background: The Z/^* — > r+r 
background is estimated using the r-embedded Z/j* — t- sample outlined in sec- 

tion 6.1. The jet activity in the embedded events is independent of the Z boson decay 
mode. Taking advantage of this feature, the embedding sample is also used to validate 
the simulated Z/^* — )• e^e~ and — )• background samples for the correct 6-jet 

fraction, which may affect the background estimation after imposing the 6-tag requirement. 
Correction factors are derived by comparing r-embedded Z/^* — t- ^'^ ^~ events with simu- 
lated Z/^* —7- T+r" events before and after the 6-tagged sample selection. The correction 
factors are calculated to be /^^ = 1.08 it 0.23 and /^^ = 1.11 it 0.13 for the electron and 
muon channels, respectively, where the quoted uncertainties are statistical. The effect on 
these correction factors from the ti contribution in the control region is studied by remov- 
ing the 50 GeV maximum Et requirement on the 6-jet. A 7% systematic uncertainty is 
obtained. These factors are applied to the simulated Z/^* — )• e^e~ and Z/'^* — >• fi^n" 
background samples passing the 6-tagged sample selection. 

Estimation of the tt background: The simulated ti samples are normalised from data 
using a top-enriched control region. This control region is defined by applying the rjcpThad 
selection criteria up to the requirement of a erhad or ^Thad pah in the event, with no 
requirement on the transverse mass. The highest-£'T jet in the event must be identified as 
a 6-jet, with Et in the range 50 GeV to 150 GeV, and a second highest-E'T jet must satisfy 
the same 6-jet identification requirement. This results in a control region with a purity of 
tt events over 90% and negligible signal contribution. The ti correction factor is derived 
in a manner similar to that of the W + jets correction factor, and a value of /^j = 0.88 
±0.04 (stat.) ±0.14 (syst.) is obtained with the systematic uncertainty due primarily to 
the 6-jet identification efficiency. 

Estimation of the multi-jet background: For the multi-jet background estimation, 
the ABCD method is used by defining four regions according to whether the charge of the 
r jet and lepton have opposite sign or same sign, and whether the selected lepton passes 
or fails the isolation criteria. These requirements are summarised in table 2. In regions C 
and D the contribution from processes other than the multi-jet background is negligible, 
while in region B there is a significant contribution from other backgrounds, in particular 
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Figure 4. MMC mass distributions for tlic h/A/H — > TiopThad final state. The MMC mass, 
j^MMC^ is shown for the 6-tagged (Icft-liand side) and 6- vetoed samples (right-hand side) for the 
combined TeThad and r^Thad samples. The data are compared to the background expectation and a 
hypothetical MSSM signal (itla = 150 GeV and tan/3 = 20). The background uncertainties include 
statistical and systematic uncertainties. The contributions of the diboson and VF+jets background 
processes are combined and labelled "Other electroweak" . 



Z/7* -|-jets and W + jets, which is subtracted from the data sample using estimates from 
simulation. The systematic uncertainty on the predicted event yield is estimated by varying 
the definitions of the regions used, and by testing the stability of the r^/^) ratio across the 
^MMC j-g^j^ge. The resulting uncertainty is 7.5% in the r^Thad channel and 15% in the 
TeThad channel. 

Results: The number of observed TicpT^ad events in data, along with predicted event 
yields from background processes, are shown in table 4. The observed event yields are 
compatible with the expected yields from Standard Model processes within the uncertain- 
ties. The MMC mass distributions for these events, with TeThad and TgThad statistically 
combined, are shown in figure 4. 

6.4 The h/A/H —> Thad'Hiad decay channel 

Signal topology and event selection: Events in this channel are selected by a di-Thad 
trigger with transverse momentum thresholds of 29 GeV and 20 GeV for the two Thad can- 
didates. Events containing identified electrons or muons with transverse momenta above 
15 GeV or 10 GeV, respectively, are vetoed. These vetoes suppress background events and 
ensure that the channels are statistically independent. Two Thad candidates with opposite- 
sign charges are required, one passing the tight Thad identification requirements and the 
second passing the medium criteria. These two leading Thad candidates are required to 
match the reconstructed Thad trigger objects each within a cone with radius parameter 
AR < 0.2. These two Thad candidates are required to have transverse momenta above 
45 GeV and 30 GeV, respectively. These values are chosen such that the plateau of the 
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Muon Channel (r^Thad) 
6-tagged sample 6-vetoed sample 





Z/7* ^ r+T" 




86 


±15 


4800 


±700 




+ jets 




19 


+6 

-8 


780 


+100 
-140 
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+5 
-4 


350 


+100 
-90 




Top 




14.5 


+3.5 
-2.7 
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+20 
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Diboson 




0.8 


±0.4 


38 


+6 
-5 




Multi-Jet 




51 


±11 
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+140 
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Total 




180 


±20 


6600 


±800 




^nal mA = 150GeV,tan/3 = 
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bb{h/A/H^ tt) 




20 


+5 

-6 
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+27 
-35 




gg ^h/A/H^ tt 




1.2 


±0.6 


115 


±16 




Data 
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6-tag 


;ged sample 


6-vetoed sample 




Z/7* ^ r+T- 




40 


±20 


2700 


±500 




+ jets 




18 


+9 

-12 


740 


+110 
-160 








19 


±10 


700 


+350 
-270 




Top 




15 


±3 


106 


+20 
-21 




Diboson 




1.0 


+0.4 
-0.5 


29 


+5 
-4 




Multi-Jet 




60 


±15 


920 


+230 
-240 




Total 




150 


±30 


5200 


±600 




i;nal tua = 150GeV,tan/3 = 


20 












hb{h/A/H^ tt) 




15 


+3 
-5 


138 


+22 
-29 




gg ^h/A/H^ tt 




1.2 


+0.6 
-0.4 


99 


+15 
-14 




Data 




175 




5034 





Table 4. The number of events observed in data and the expected number of signal and back- 
ground events for the h/A/H — >■ riopThad channel. Simulated event yields are normalised to the 
integrated luminosity of the data sample, 4.7 fb~^. The predicted signal event yields correspond 
to the parameter choice niA = 150 GeV and tan /3 = 20 and include both the 6-associated and the 
gluon-fusion production processes. 
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trigger turn-on curve is reached and the electroweak and multi-jet backgrounds are sup- 
pressed effectively. The missing transverse momentum is required to be above 25 GeV 
to account for the presence of neutrinos originating from the r decays and to suppress 
multi-jet background. 

The selected events are split into a 6-tagged sample and a 6-vetoed sample to exploit 
the two dominant production mechanisms for neutral Higgs bosons in the MSSM. Events 
in which the leading jet is identified as a b-jet are included in the 6-tagged sample. The 
transverse momentum of this jet is restricted to the range of 20 GeV to 50 GeV to reduce 
the tt background. Events without jets, or in which the leading jet is not identified as 
a 5-jet, are included in the 6-vetoed sample. Due to the higher background levels in this 
sample, the threshold on the transverse momentum of the leading Thad candidate is raised 
to 60 GeV. 

Identification efficiency and misidentification corrections for hadronic r de- 
cays: The Thad identification efficiencies, the Thad trigger efficiencies and the correspond- 
ing misidentification probabilities are corrected for differences observed between data and 
simulation. For the di-Thad trigger it is assumed that these identification and misidenti- 
fication efficiencies can be factorised into the efficiencies of the corresponding single- Thad 
triggers with appropriate transverse momentum requirements. This factorisation is vali- 
dated using a simulated event sample. The single- Thad trigger efficiency for real hadronically 
decaying r leptons with respect to the offline Thad selection was measured using a tag-and- 
probe analysis with Z — )• r^Thad data. A correction factor for the simulation was derived 
as a function of the transverse momentum of each of the two Thad candidates, and each 
event was weighted by the product of these factors. The probability to misidentify a jet as 
3- Thad is extracted for both the trigger and the Thad identification algorithm by analysing 
jets in a high-purity W{-^ fxu)+^ets sample. A correction factor derived on the basis of 
these probabilities is applied to the simulation when a jet is misidentified as a Thad- The 
statistical and systematic uncertainty of these correction factors leads to an uncertainty of 
21% on the W + jets background and 4% to 6% on the tt background. 

Estimation of the Zf-y* t~^t~ and W{—^ T^')-|-jets backgrounds: The esti- 
mates of the Z/^* —7- r~'"r~ and VF(— > rzv)+jets backgrounds are taken from simulation 
and are validated using r-embedded Z/^* — )• and W{-^ ;Ui/)+jets samples. The r- 

embedded Z/^* — ?• fJ-~^fJ,~ data are used to validate the simulation, rather than to provide 
the main estimate for the Z/j* — t- t'^t~ background, in the h/A/H — t- Thad^had chan- 
nel. The di-Thad trigger is not modelled in the r-embedded Z/y* — t- t^t~ data, making 
it difficult to apply the embedding technique. Correction factors for the efficiency of the 
6-jet identification requirement on the leading jet are derived in a way equivalent to that 
described in section 6.3. For the Z/^* +jets background a factor of fzb = 1.24ib0.34 (stat.) 
is derived by comparing the simulated and embedded Z/y* — )■ t^t~ samples. With the 
embedded W{-^ Tz/)+jets sample, no such correction factor may be derived in this way as 
the contamination from ti events is quite significant once the 6-jet identification require- 
ments are applied. Instead, the procedure in section 6.3 is applied with the selection of 
this channel. A correction factor of fy/b = 1.00 ±0.31 (stat.) is derived for W{-^ rz/)+jets 
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Region Charge correlation Hadronic r decay identification requirement 



A (Signal Region) 



Opposite sign 
Same sign 

Opposite sign 
Same sign 



Pass 
Pass 
Fail 
Fail 



B 
C 
D 



Table 5. Control regions for the estimation of the multi-jot background for the h/A/H Thad'Hiad 
selection: Events are categorised according to the product of the electric charges of the two Thad 
candidates and the Thad identification. "Pass" refers to one Thad candidate passing the tight and 
the other Thad candidate passing the medium identification selection. "Fail" refers to all events in 
which the two Thad pass the loose identification selection but do not satisfy the selection of the 
"Pass" category. 

events. The correction factors for the 6- vetoed sample are found to be close to 1 and 
uncertainties are negligible, hence no correction is applied. 

Estimation of the multi-jet background: The multi-jet background is estimated us- 
ing the ABCD method, by splitting the event sample into four regions based on the charge 
product of the two leading Thad candidates and whether the nominal Thad identification 
requirements of these two Thad candidates are met. These variables can be assumed to be 
uncorrelated for multi-jet events. Table 5 illustrates the definition of the four regions. 

The shape of the MMC mass distribution for the multi-jet background in the signal 
region is taken from region C for the 6-vetoed sample and from region B for the 6-tagged 
sample. Contributions from electroweak backgrounds are subtracted from each control 
region using simulation. The uncertainties on these backgrounds lead to uncertainties on 
the multi-jet estimate of 7% for the 6-tagged sample and 5% for the 6-vetoed sample. 

Results: The number of observed ThadT'had events in data, along with predicted event 
yields from background processes, is shown in table 6. The observed event yields are com- 
patible with the expected yields from Standard Model processes within the uncertainties. 
The MMC mass distributions for these events are shown in figure 5. 

7 Systematic uncertainties 

Data-driven background estimation: Where possible, event yields and mass distri- 
butions for the background are estimated using control samples in data. The specific 
techniques and their associated uncertainties have been presented in the relevant sections. 
The effect of these uncertainties on the predicted background event yield is less than 5% for 
the jJL^ lJi~ channels and less than 15% for the t^t~ channels, and is usually small compared 
to the systematic uncertainties from the simulated samples. 

Cross-section for signal and background samples: The uncertainties on the signal 
cross-sections are estimated to be 10% to 20%, depending on the values of and tan/3, 
for both gluon- fusion and 6-associated Higgs boson production [29]. An uncertainty of 5% 
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Table 6. The observed number of events in data and the expected number of signal and background 
events for the h/A/H — )■ Thad'Hiad channel. Simulated event yields are normalised to the total 
integrated luminosity, 4.7 fb~^. The data arc compared to the background expectation and an 
added hypothetical MSSM signal (m^ = 150 GeV and tan/3 = 20). Because of the subtraction of 
the Z /^* — T~^T~ background in the control regions used in the multi-jet background estimation 
the uncertainties of the Z/^* — > r+r^ and the multi-jet background are anti-correlatcd. 



350 



C5 



0) 

> 

LU 



, b-tagged sample 




• Data 2011 

m„=1 50 QeV, tan p=20 

I I Multi-jet 

Z — > TT 
P~1 W^TV 

Top 

^ Bkg. uncertainty 



b-vetoed sample 



100 200 300 400 500 600 




• Data 2011 

mj=150 GeV, tan p=20 

I I IVIuiti-jet 
EllZ^Tx 

I I W^TV 

Bkg. uncertainty 
ATLAS 

Ji = 7TeV,|Ldt = 4.7fb"' I 



200 400 600 800 1000 1200 



nMMC 



[GeV] 



[GeV] 



Figure 5. MMC mass distributions for the h/A/H — >• rhad''"had final state. The MMC 
mass , m/f/^^ , is shown for the ^-tagged (left-hand side) and 6- vetoed samples (right-hand side). 
The data are compared to the background expectation and an added hypothetical MSSM signal 
[mA = 150 GeV and tan/3 = 20). The background uncertainties include statistical and systematic 
uncertainties. 



is assumed on the cross-sections for W and Z boson background production [63, 64]. Uncer- 
tainties due to the parton distribution functions and the renormalisation and factorisation 
scales are included in these estimates. 
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Acceptance modelling for simulated samples: The uncertainty on the acceptance 
from the parameters used in the event generation of signal and baclcground samples is also 
considered. This is done by evaluating the change in kinematic acceptance after varying the 
relevant scale parameters, parton distribution function choices, and if applicable, conditions 
for the matching of the partons used in the fixed order calculation and the parton shower. 
Furthermore, the effects of different tunes of the underlying event activity are considered. 
The resulting uncertainties are typically 2% to 20%, depending on the sample and the 
channel considered. 

Electron and muon identification and trigger: The uncertainties on electron or 
muon trigger and identification efficiencies are determined from data using samples of W 
and Z decays [50, 51]. The trigger efficiency uncertainties are below about 1%. The 
identification efficiency uncertainties are between 3% and 6% for electrons and below 1.8% 
for muons. The total effect of these uncertainties on the event yield is no greater than 4% 
in any channel. 

Hadronic r identification and trigger: The uncertainties related to hadronic r trigger 
and identification efficiencies are also studied with data [58] . The di-Thad trigger efficiency 
used in the rhad''"had channel has an uncertainty of 2% to 7%. The identification efficiency 
uncertainty for hadronic r decays is about 4% for reconstructed Thad Pt above 22 GeV and 
8% otherwise. These uncertainties are most important in the Thad'^had channel, where the 
effect on the estimated signal yield reaches 11%. 

6-jet identification: The 6-jet identification efficiency and the misidentification prob- 
abilities for jets other than 6-jets have been measured in data [57, 65]. The associated 
uncertainties are treated separately for all jet flavours and depend on jet Et and r]. Typ- 
ical uncertainty values are around 5% for 6-jets and between 20% and 30% for other jets, 
leading to a total uncertainty close to 5% on the event yield in the 6-tagged channels. 

Energy scale and resolution: The uncertainty on the acceptance due to the energy 
measurement uncertainty in the calorimeter is considered for each identified object corre- 
sponding to the clusters in the calorimeters. For the clusters identified as electrons, typ- 
ically a 1% (3%) energy scale uncertainty is assigned for the barrel (end-cap) region [50]. 
The energy scale uncertainties for clusters identified as hadronic r decays and jets are 
treated as being fully correlated, and are typically around 3% [56, 58]. The uncertainty 
in the muon energy scale is below 1%. The acceptance uncertainty due to the jet energy 
resolution, which affects the pT thresholds used to define the 6-tagged and 6- vetoed sam- 
ples, is typically less than 1%. The systematic uncertainty due to the energy scales of 
electrons, muons, hadronic r decays and jets is propagated to the ii^™'^^ vector. Additional 
uncertainties due to different pile-up conditions in data and simulation are also considered. 

The uncertainty on the acceptance due to the energy scale and resolution variations 
reaches up to 37% for signal in the r^ad'^'had 6-tagged channel, but is usually less than 10% 
for channels with fewer Thad- 
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Luminosity: The simulated sample event yield is normalised to the integrated luminosity 
of the data, which is measured [66, 67] to be 4.7 fb~^ and 4.8 fb~^ for the r+r" and fj,~^n~ 
channels, respectively, and has an uncertainty of 3.9%. This is applicable to all signal and 
background processes which are not normalised using data-driven methods. 



8 Statistical analysis 

The statistical analysis of the data employs a binned likelihood function. Each one of the 
///i, TeT^, TeThad) TfiTha.d and Thad''"had final statcs is split iuto a 6-tagged and a 6-vetoed 
sample. The likelihood in each category is a product over bins in the distributions of the 
MMC mass in the signal and control regions. 

The expected number of events for signal (sj) and background (bj), as well as the ob- 
served number of events (Nj) in each bin of the mass distributions, enter in the definition 
of the likelihood function C{^,9). A "signal strength" parameter (/x) scales the expected 
signal in each bin. The value /.i = corresponds to the background-only hypothesis, while 
/i = 1 corresponds to the signal-plus-background hypothesis with all Higgs bosons having 
the masses and cross-sections specified by the point considered in the m^-tan /3 plane for 
the MSSM exclusion limit. Signal and background predictions depend on systematic un- 
certainties that are parameterised by nuisance parameters, 6, which in turn are constrained 
using Gaussian functions, J-'q, so that 

£ {fi, e)= W Tp {Nj \^.sj+bj)ll Tg {9i I 0, 1) , (8.1) 

j = bin and 8i 
category 

where {Nj \ ■ Sj + bj ) denotes the Poisson distribution with mean fj, ■ Sj + bj for vari- 
able Nj. The correlations of the systematic uncertainties across categories are taken into 
account. The expected signal and background event counts in each bin are functions of 
0. The parameterisation is chosen such that the rates in each channel are log- normally 
distributed for a normally distributed 6. 

To calculate the upper limit on /_i for a given signal hypothesis, the compatibility of 
the observed or expected dataset with the signal-plus-background prediction is checked 
following the modified frequentist method known as CL^ [68]. The test statistic q^, used 
in the upper limit derivation, is defined as 



9m 



if fi > fi, 



where fi and refer to the global maximum of the likelihood and 0^ corresponds to the 
conditional maximum likelihood for a given The asymptotic approximation [69] is used 
to evaluate the probability density functions rather than performing pseudo-experiments; 
the procedure has been validated using ensemble tests. 

The significance of an excess in data is quantified with the local pQ-value, the proba- 
bility that the background processes can produce a fiuctuation greater than or equal to the 
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excess observed in data. The test statistic qo, which is used in the local po- value calculation, 
is defined as: 

r_21n(^) ifA>0, 
qo=l ^ ^('^''') ^ (8.3) 

if/i<0, 

where the notation is the same as in Equation 8.2. The equivalent formulation in terms 
of the number of standard deviations (a), Zq, is referred to as the local significance and is 
defined as: 

Zo = $-i(l-po), (8.4) 

where is the inverse of the cumulative distribution of the Gaussian distribution. The 
local pQ-value is estimated using the asymptotic approximation for the go distribution [69] . 



9 Results 

No significant excess of events above the background-only expectation is observed in the 
considered channels in 4.7 fb~^ to 4.8 fb~^ of ^/s = 7TeV proton-proton collision data. A 
95% CL upper limit on tan f3 is set for each point using the frequentist method described 
in section 8. This is done using Higgs boson cross-sections calculated in the m^^^ scenario 
with n > [11]. Results for each of the jifi, TeT^, ricpThad and Thad''"had final states, as well 
as for their statistical combination, can be seen in figure 6. The tightest constraint is at 
ruA = 130 GeV, where values of tan/3 > 9.3 are excluded. The expected exclusion for the 
same point is tan /? > 10.3. The exclusion of parameter space is significantly increased in 
comparison to earlier results by the ATLAS Collaboration [7] and complementary to the 
excluded parameter space from searches at LEP [13]. A significant portion of the MSSM 
parameter space that is not excluded is still compatible with the assumption that the newly 
discovered particle at the LHC is one of the neutral CP-even MSSM Higgs bosons [19, 20]. 
The lowest local po-values per channel are 0.014 {2.2a) at ruA = 125 GeV for the fi^ channel, 
0.014 (2.2(7) at ruA = 90 GeV for the TeT/, channel, 0.067 (1.5o-) at tua = 90 GeV for the 
''"icp'^'had channel and 0.097 (1.3o") at niA = 140 GeV for the Thad^had channel. The lowest 
local pq for the statistical combination of all channels is 0.004 (2.7o") at uia = 90 GeV. The 
significance of this excess is below 2a, after considering the look-elsewhere effect in the 
range 90 GeV < niA < 500 GeV and 5 < tan/3 < 60 [70]. 

The outcome of the search is further interpreted in the generic case of a single scalar 
boson (p produced in either the gluon-fusion or 6-associated production mode and decaying 
to fJ-'^fM~ or T"^r~. Figure 7 shows 95% CL limits based on this interpretation. The 
exclusion limits for the production cross-section times the branching ratio for a Higgs 
boson decaying to or t~^t~ are shown as a function of the Higgs boson mass. 



10 Summary 

A search is presented for the neutral Higgs bosons of the Minimal Supersymmetric Stan- 
dard Model in proton-proton collisions at a centre-of-mass energy of 7 TeV with the ATLAS 
experiment at the LHC. A significant portion of the available MSSM parameter space is 
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Figure 6. Expected (dashed line) and observed (solid line) 95% CL limits on tan^ as a function 
of uiA for the statistical combination of all channels along with the ilcr and ±2 a bands for the 
expected limit are shown on the left plot. Values of tan/3 greater than the shown lines are excluded. 
The 95% CL limits for the expected limit (dashed lines) and the observed limit (continuous lines) for 
each of the fxfi, TeT^, riopThad and Tiiad''"had channels and their statistical combination are shown on 
the right plot. The 95% CL exclusion region from neutral MSSM Higgs boson searches performed 
at LEP[13] is shown in a hatched style. 
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Figure 7. Expected (dashed line) and observed (solid line) 95% CL limits on the cross-section for 
gluon-fusion and 6-associated Higgs boson production times the branching ratio into r and fi pairs, 
respectively, along with the ±la and ±2 a bands for the expected limit. The combinations of all 
rr and fifi channels are shown. The difference in the exclusion limits obtained for the gluon-fusion 
and the 5- associated production modes is due to the different sensitivity from the 6-tagged samples. 

consistent with the assumption that the newly discovered particle at the LHC is one of 
the neutral CP-even MSSM Higgs bosons. The study is based on a data sample that 
corresponds to an integrated luminosity of 4.7 fb~^ to 4.8 fb~^. The decay modes of the 
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Higgs bosons considered are h/A/H — )• j-i^fi^, h/A/H — )• TeT^ , h/A/H — )• riepThad and 
h/A/H — )• Thad''"had- The analysis selection criteria exploit the two main production mech- 
anisms in the MSSM, the gluon-fusion and 6-associated production modes, by introducing 
categories for event samples with and without an identified 6-jet. Since no excess of events 
over the expected background is observed in the considered channels, 95% CL limits are 
set in the m^-tan j3 plane, excluding a significant fraction of the MSSM parameter space. 
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G. Gonzalez Parra^^, M.L. Gonzalez Silva^"^, S. Gonzalez-Sevilla^^, J.J. Goodson^^^, 
L. Goossens^°, P.A. Gorbounov^^, H.A. Gordon^^, I. Gorelov^^^, G. Gorfine^'^^, 

B. Gorini^o, E. Gorini'^2a,72b^ Gorisek'^^ E. Gornicki^s, A.T. Gosliaw^ M. Gosselink^os^ 
M.I. Gostkin^^, I. Gough Eschrich^^s ^ Gouighri^^Sa^ j3_ Goujdamii^Sc^ M.P. Goulette^^ 
A.G. Goussiou^^s, C. Goy^ S. Gozpinar^s, I. Grabowska-Bold^^ , P. Grafstrom20^'20b ^ 
K-J. Grahn^^, E. Gramstad^^^, F. Grancagnolo^^** , S. Grancagnolo^^ , V. Grassi^^^, 

V. Gratchev^^i^ jj.M. Gray^o, J.A. Grayl'^^ E. Graziani^^^a^ q.G. Grebenyuk^^i^ 

T. Greenshaw'^'^ , Z.D. Greenwood^^'", K. Gregersen^^, I.M. Gregor^^, P. Grenier^^^, 
J. Griffitlls^ N. Grigalashvili^^, A.A. Grillo^^?^ Grimm^\ S. Grinstein^^ Ph. Gris^^, 
Y.V. Grishkevich^^ J.-F. Grivazl^^ A. Grohsj ean^^, E. Gross^'''^, J. Grosse-Knetter^^, 
J. Groth-Jensen^^^ K. Grybel^^\ D. Guestl■^^ C. Guicheney^^, E. GuidoSO^'^o*^, 
T. Guilleminiis^ S. Guindon^^, U. GuF^ J. Gunther^^s, B. Guo^^^^ J. Guo^^, 
P. Gutierrez^", N. Guttmani^s, Q. Gutzwiller^'^^^ C. Guyot^^e^ q Gwenlan^^^^ 

C. B. Gwilliam'^^ A. Haasl°^ S. Haas^o, C. Haber^^ H.K. Hadavand^, D.R. Hadley^^ 
P. Haefner^i, F. Hahn^o, Z. Hajduk^^, H. Hakobyani'^^ D. Halll^^ G. Halladjian62, 
K. Hamacher^'^^, P. Hamal^^^, K. Hamano^^, M. Hamer^^, A. Hamilton^''^'^'^, 

S. Hamiltoni6\ L. Han^^b, K. Hanagaki^^^ K. Hanawa^^^, M. Hancel^ C. HandelS\ 

P. Hanke^^'^, J.R. Hansen^^, J.B. Hansen^^, J.D. Hansen^^, P.H. Hansen^^, P. Hansson^^^, 

K. Hara^^'^, T. Harenberg^^^ , S. Harkuslia^'^, D. Harper^^, R.D. Harrington^^ , 

0. M. Harris^^^, J. Harterf*^, F. Hartjes^°^, T. Haruyama'^^ , A. Harvey^*^, S. Hasegawa^°\ 
Y. Hasegawa^^°, S. Hassani^^e^ Haug^^, M. Hauschild^" , R. Hauser*^^, M. Havranek^^ 
CM. Hawkes^*^, R.J. Hawkings^o, A.D. Hawkins^^ T. Hayakawa^^, T. Hayaslii^'^o ^ 

D. Hayden^s^ CP. Hays^l^ H.S. Hayward^^^ S.J. Haywood^^g^ g Head^^^ V. Hedberg^^ 
L. Heelan^, S. Heim^^^, B. Heinemann^^ , S. Heisterkamp^^, L. Helary^^, C. Heller^^, 

M. Heller^o, S. Hellmani'^Sa.Meb^ Hellmicli^i, C. Helsens^^^ R.C.W. Henderson7\ 
M. Henke^'^'^, A. Henrichs^''^, A.M. Henriques Correia^°, S. Henrot-Versille^^^, 
C. Hensel^^, CM. Hernandez^, Y. Hernandez Jimenez^^^, R. Herrberg^^, G. Herten^^, 
R. Hertenberger^^, L. Hervas^°, G.G. Hesketh^'^, N.P. Hessey^°^, R. Hickling^^, 

E. Higon-Rodriguez^^^ J.C HilF, K.H. Hiller^^ S. Hillert^^ S.J. Hillieri^ 

1. Hinchliffe^^ E. Hines^^o^ Hirosel^^ F. Hi^sch'^^ D. Hirscllbuehll^^ J. Hobbs^''^ 
N. Hodl5^ M.C Hodgkinson^^g^ p_ Hodgson^^s, A. Hoecker^o, M.R. Hoeferkampi°^ 
J. Hoffman^o, D. Hoffmann^^^ M. HohlfeldS\ M. Holderi4\ S.O. Holmgren^^^a^ 

T. Holyi26, J.L. Holzbauer^s, T.M. Hong^^o, L. Hooft van Huysduynen^os, S. Horner4^ 
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J-Y. Hostachy^^, S. Hou^^\ A. Hoummada^^^'', J. Howard^^'^, J. Howarth*^^, I. Hristova^^, 
J. Hrivnac^l^ T. Hryn'ova^, P.J. Rsu^\ S.-C. Hsu^^s^ ^ jju^^, Z. Hubacek^o, 

F. Hubaut^^, F. Huegging^^ A. Huettmann^^ ^ Huffman^^^, E.W. Hughes^^, 

G. Hughes^\ M. Huhtinen^o, M. Hurwitz^^, N. Huseynov^^''', J. Huston^^, J. Huth^'^, 

G. lacobucci^^, G. lakovidis^'^, M. Ibbotson®^, I. Ibragimov^'^^ , L. Iconomidou-Fayard^^^, 
J. Idarraga^^^, P. lengo^"^^, O. Igonkina^^^, Y. Ikegami^^, M. Ikeno^^, D. Iliadis^^^, 
N. Ilicl^^ T. Ince^^ P. loannou^ M. Iodicei34a^ lordanidou^ V. Ippolito^32a,i32b^ 
A. Irles Quiles^'^'^, C. Isaksson^'^^, M. Ishino'^'^, M. Ishitsuka^^'^, R. Ishmukhametov^°'^ , 

C. Issever^l^ S. Istin^^^, A.V. Ivashin^^s^ Iwanski^^, H. Iwasaki^s, J.M. Izen^\ 
V. Izzo^°2a^ g Jackson^^o^ Jackson■^^ P. Jackson^ M.R. JaekeF, V. Jain^, 

K. Jakobs''^ S. Jakobsen^^, T. Jakoubek^^s, j. Jakubek^^e^ q. Jamin^si^ Jana^l^ 
E. Jansen^^, H. Jansen^'', J. Janssen^^, A. Jantsch^^, M. Janus^*^, R.C. Jared^^^, 
G. Jarlskog'''^ , L. Jeanty^^, I. Jen-La Plante'^^, G.-Y. Jeng^^°, D. Jennens^^, P. Jenni^°, 
A.E. Loevschall-Jensen36, P. Jez^e, S. Jezequel^ M.K. Jha^O'^, H. Ji^^^ W. Ji^^, J. Jial^^ 
Y. Jiang^^*^, M. Jimenez Belenguer^^, S. Jin^^*^, O. Jinnouchi^^^, M.D. Joergensen^^ , 

D. Joffe'^°, M. Johansen^^^^'^^^'^, K.E. Johansson^^^a^ p johansson^^g^ Johnert''^ 
K.A. Johns^ K. Jon-And^'^'^'^'i'''^'', G. Jones^^^, R.W.L. Jones^\ T.J. Jones^^ C. Joram^o, 
P.M. Jorge^24a^ j3 joshi82^ j_ jovicevic^^^ T. Jovin^^b^ x. Ju^■^^ C.A. Jung^^ 

R.M. Jungst^°, V. Juranek^25^ p_ juggelGi^ j^g^^ Rozas^^, S. Kabana^'^, M. Kaci^'^'^, 
A. Kaczmarska^^, P. Kadlecik^^, M. Kado^^^, H. Kagan^^^, M. Kagan^^, 

E. Kajomovitz^^^, S. Kalinin^''^, L.V. Kalinovskaya^^ , S. Kama""^, N. Kanaya-^^^, 

M. Kaneda^o, S. Kaneti^^, T. Kanno^^'^, V.A. Kantserov^^, J. Kanzaki^^, B. Kaplan^o^, 

A. Kapliy^^, D. Kar^^, M. Karagounis^^ , K. Karakostas^", M. Karnevskiy^^'^ , 

V. Kartvelishvili'^\ A.N. Karyukhin^^s, L. Kashif^^^ G. Kasieczka^^'^ , R.D. Kass^^^ 

A. Kastanas^^, M. Kataoka^, Y. Kataoka^^^, J. Katzy''^, V. KaushiF, K. Kawagoe'^^, 

T. Kawamoto^^^ , G. Kawamura^^, S. Kazama^^^, V.F. Kazanin^"^, M.Y. Kazarinov^^, 

R. Keeler^SQ^ p-p^ Keener^^o^ ^ Kehoe''°, M. Keil^^, G.D. Kekelidze^^, J.S. Keller^^s^ 

M. Kenyon^'^, H. Keoshkerian^ , O. Kepka^^^, N. Kerschen'^^, B.P. Kersevan'''^ , 

S. Kersten^■^^ K. Kessoku^5^ J. Keung^^^ F. Khalil-zadai\ H. Khandanyan^^^'^'i^'''^, 

A. Khanov"^ D. Kharchenko''^ A. Khodinov^^ A. Khomich^^'^, T.J. Khoo^^, 

G. Khoriauli^^, A. Khoroshilov^^^, V. Khovanskiy^^ , E. Khramov^^, J. Khubua^^*^, 

H. Kim^46a,i46b^ g jj_ Kimi^o^ N. Kimura^^i, O. Kmdl^ B.T. King^^ M. King^e, 
R.S.B. Kingii^, J. Kirk^^g^ Kiryunin^s, T. Kishimoto^^ , D. Kisielewska^^ , 

T. Kitamura'56^ T. Kittelmann^^s, K. Kiuchii^o, E. Kladiva^'^^'^ , M. Klein^^ U. Klein^^ 

K. Kleinknecht^i, M. Klemetti8^ A. Klier^^^^ P. Klimek^^ea.ueb^ Klimentov^s, 

R. Klingenberg''^, J. A. Klinger^^ E.B. Klinkby3^ T. Klioutchnikova^^ , P.F. Klok^°^, 

S. Klousl°^ E.-E. Kluge^'^'^, T. Kluge■^^ P. Kluitl°^ S. Kluth^^ E. Kneringer6\ 

E.B.F.G. Knoops^^^ A. Knue^^, B.R. Ko''^ T. Kobayashi^5^ M. Kobel'^^ M. Kocian^^^^ 

P. Kodys^^T^ ^ Koneke^", A.C. Konig^O'^, S. Koemg'^\ L. K6pkeS\ F. Koetsveld^°^, 

P. Koevesarki^i, T. Koffas^^, E. Koffeman^os ^ p.A. Kogan^l^ S. Kohlmann^7^ F. Kohn^^, 

Z. Kohout^26^ rp_ Kohriki^s, T. Koi^^^^ G.M. Kolachev^o^'*, H. Kolanoski^^^ 

V. Kolesnikov^^, I. Koletsou^^'^, J. KolF, A.A. Komar9^ Y. Komorii^s, T. Kondo'^^ 

T. Kono''^''', A.I. Kononov^^, R. Konoplich^°^'*, N. Konstantinidis'^'^, R. Kopeliansky^^^^ 
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S. Koperny^s, K. KorcyF, K. Kordas^^'^, A. Kom^l^ A. Korol^o^ I. Korolkov^^ 
E.V. Korolkova^^^ , V.A. Korotkov^^®, O. Kortner^^, S. Kortner^^, V.V. Kostyukhin^^, 
S. Kotov^s, V.M. Kotov'5^ A. Kotwal^^^ C. Kourkoumelis^ V. Kouskoura^5^ 
A. Koutsman^ssa^ ^ Kowalewski^^g ^ T ^. Kowalski^^, W. Kozanecki^^e , A.S. Kozhini^s, 

V. Kral^^^, V.A. Kramarenko^^, G. Kramberger''^, M.W. Krasny'^^, A. Krasznahorkay^'^^ , 
J.K. Kraus^^, A. Kravchenko^^ , S. Kreiss^°^, F. Krejci^^^, J. Kretzschmar'^^ , 
K. Kreutzfeldt^^, N. Krieger^"^, P. Krieger^^®, K. Kroeninger^^, H. Kroha^^, J. Kroll^^°, 
J. Kroseberg^^, J. Krstic^^'^, U. Kruchonak^^, H. Kriiger^^, T. Kruker^^, N. Krumnack^^, 
Z.V. Krumshteyn6^ M.K. Kruse^^ T. Kubota^^ S. Kuday'^'^, S. Kuehn''^ A. Kugel^*^ 
T. Kuhl^^ V. Kukhtin^^, Y. Kulchitsky'^o, S. Kuleshov32b, M. Kuna■^^ J. Kunkle^^o^ 
A. Kupcoi25, H. Kurashige^e, M. Kurata^^^, Y.A. Kurochkin^o, V. Kus^^s, 

E. S. Kuwertzi'^^ M. Kuze^^^ J. Kvitai'^^ R. Kwee^^, A. La Rosa''9, L. La Rotonda^^^'^^b^ 
L. Labarga^o, S. Lablak^^s^^ c. Lacasta^^^ F. Lacavai32a,i32b^ j_ Lacey^^, H. Lacker^^^ 

D. Lacour"^*^, V.R. Lacuesta^^^, E. Ladygin^^, R. Lafaye^, B. Laforge^^, T. Lagouri^'^^, 
S. Lai^^, E. Laisne^^, L. Lambourne^^, C.L. Lampen^, W. Lampl"^, E. Lancon^^^, 
U. Landgraf^^ M.P.J. Landon^^ V.S. Lang^^^, C. Lange^^^ A.J. Lankford^^^ F. Lanni^^, 
K. Lantzsch3°, A. Lanza^^^'', S. Laplace^^, C. Lapoire^^, J.F. Laporte^^*^, T. Lari^^'", 

A. Larner^^*^, M. Lassnig^^, P. Laurelli^'^, V. Lavorini^'^'^'^'^^, W. Lavrijsen-'^^, P. Laycock'''^, 
O. Le Dortz^*^, E. Le Guirriec^^, E. Le Menedeu^^, T. LeCompte'^, F. Ledroit-GuilW^, 
H. Leel°^ J.S.H. Lee^l^ S.C. Lee^^^, L. Lee^^^ M. Lefebvre^6^ M. Legendre^^e^ 

F. Legger^^, C. Leggett^^, M. Lehmacher^-^ , G. Lehmann Miotto'^^, A.G. Leister-'^'^^, 
M.A.L. Leite^^'^, R. Leitner^^'', D. Lellouch^'^^, B. Lemmer^'', V. Lender mann^^'^, 
K.J.C. Leney^'^5'^, T. Lenz^os^ q Lenzen^^s^ B. Lenzi^o, K. Leonhardf*^, S. Leontsinisi°, 
F. Lepold^S'^, C. Leroy93^ J-R. Lessard^^^ C.G. Lester^s, CM. Lester^^o^ j_ Leveque^ 
D. Levin^^, L.J. Levinson^'^^ ^ Lewis^^^, G.H. Lewis^°^, A.M. Leyko^^, M. Leyton^'^^ 

B. Li33b^ B. U^^, H. Lil4^ H.L. Li^^, S. Li^^b.", x. LF, Z. Liang^^^''', H. Uao^^, 
B. Liberti^33a^ p_ Lichard^^, K. Lie^^^ , W. Liebigl^ C. Limbach^i, A. Limosani^^ 

M. Limper''^^ S.C. Lin^^^'"', F. Lmdel°^ J.T. LmnemannS^ E. Lipeles^^o^ Lipniacka^'*, 

T.M. Uss^^^, D. Lissauer^s, A. Lister^^ A.M. Litke^^^ D. Liu^^^, J.B. Liu^^b^ l_ Liu^?^ 

M. Liu^^b, Y. Liu^^b, M. Livan^^^^'^^^'^, S.S.A. Livermore^l^ A. Lleres^^ 

J. Llorente Merino^°, S.L. Lloyd'^^, E. Lobodzinska'*^, P. Loch'^, W.S. Lockman^'^^, 

T. Loddenkoetter^i, F.K. Loebinger^^^ A. Logmov^'^6^ q Loh.^^^, T. Lohse^^, 

K. Lohwasser^^, M. Lokajicek^^^, V.P. Lombardo^, R.E. Long'^^, L. Lopes-^^^*^, 

D. Lopez Mateos^^, J. Lorenz^*^, N. Lorenzo Martinez^^^, M. Losada^^^, P. Loscutoff^^, 

F. Lo Sterzo^32a,i32b^ Lostyi^^^^'*, X. Lou^\ A. Lounis^l^ K.F. Loureiro^^^^ J. Love^, 

P.A. Love'^\ A.J. Lowe^''^'^', F. Lu^^a^ jj Lubatti^3^ C. Luci^32a,i32b^ Lucotte5^ 

D. Ludwig''^ L Ludwig^^ J. Ludwig'^^ F. Luehring^o, G. Luijckxl°^ W. Lukas6\ 

L. Luminarii32a^ ^ Lund^^^ B. Lund-Jensen^'^^ B. Lundberg^^ J. Lundbergi^^^'i^^'^, 

O. Lundberg^^6a,i46b^ Lundquist^^, M. LungwitzS\ D. Lynn^^, E. Lytken^^ H. Ma^^, 

L.L. Ma^^3, G. Maccarrone^^, A. Macchiolo^*^ , B. Macek'^'^, J. Machado Miguens^^^'', 

D. Macina^o, R. Mackeprang^^ , R.J. Madaras^^ H.J. Maddocks'^\ W.F. Mader^^, 

T. Maeno^s, P. Mattigl^^ S. Miittig^^^ L. Magnoni^^s, E. Magradze^^, K. Mahboubi^^ 

J. Mahlstedt^os^ g_ Mahmoud^^^ G. Mahoutl^ C. Maiani^^e, C. Maidantchik^^^ ^ 
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A. Maio^24a,c^ Majewski^s, Y. Makida6^ N. Makovec"^ P. Mal^^e, B. Malaescu^^ 

Pa. Malecki39, P. Malecki^^, V.P. Maleev^^i, F. Malek^s, U. Mallik^^ D. Malon^, 

C. Malone^^^, S. Maltezos^°, V. Malyshev^^'^ , S. Malyukov^^, J. Mamuzic^^'', 

A. Manabe6^ L. Mandelli^^'^, I. Mandic'^^, R. Mandrysch62, J. Maneira^24a^ 

A. Manfredini^^ , L. Manhaes de Andrade Filho^^'^, J. A. Manjarres Ramos^^^, A. Mann^^, 

P.M. Manning^^'^, A. Manousakis-Katsikakis'^ , B. Mansoulie^^*^ , R. Mantifel^^, 

A. MapellF, L. Mapelli^^, L. March^^'^, J.F. Marchand^^ F. Marchese^^^^'^^^'", 

G. Marchiori^^, M. Marcisovsky^^s ^ qj, Marino^*^'^, F. Marroquim^^^, Z. Marshall^°, 

L.F. Marti^^, S. Marti-Garcia^'^^, B. Martin^^, B. Martin^^, J.P. Martin^^, T.A. Martin^^, 

V.J. Martin^'5, B. Martin dit Latour''^, S. Martin-Haugh^^'^ , H. Martinez^^^, 

M. Martinez^^^ Martinez Outschoorn^'^, A.C. Martyniuk^'^'^ , M. Marx^^^ 

F. Marzanoi32a^ a. Marzin^^^ L. Masetti®\ T. Mashimo^^^, R. Mashinistov*^^ , 

J. Masik82, A.L. Maslennikov^o^ I. Massa20'^'20b^ g. Massarol°^ N. Massol^ 

P. Mastrandrea^^®, A. Mastroberardino^'^'^'^^^, T. Masubuchi^^^, H. Matsunaga^^^, 

T. Matsushita^'5, C. Mattravers^^^''^, J. Maurer'^^^ g Maxfield^^^ j3_A. Maximov^o^'^, 

R. Mazini^^\ M. Mazur^^ L. Mazzaferro^^^'"'^^^'^, M. Mazzanti^*^'" , J. Mc Donald^^, 

S.P. Mc Kee^^ A. McCarn^es, R.L. McCarthy^^^ T.G. McCarthy^^, N.A. McCubbin^^s, 

K.W. McFarlane^^'*, J.A. Mcfayden^^g^ q Mchedlidze^^'^ , T. Mclaughlan^^ 

S.J. McMahoni29, R.A. McPherson^^S'^, A. Meade^"^, J. Mechnichl°^ M. Mechtell^^ 

M. Medinnis^2^ g^ Meehan^^ R. Meera-Lebbai^^\ T. Meguro^^*^, S. Mehlhase^*^, 

A. Mehta^^, K. Meier^^^, B. Meirose'^^, C. Melachrinos3\ B.R. Mellado Garcia^'^^, 

F. Meloni89^'89b, l. Mendoza Navas^^^^ 2. Meng^^i.^;^ a. Mengarelli^o^'^ob, S. Menke^^ 

E. Meoni^si^ j^jvi. Mercurio^^ P. Mermod''^ L. Merola^02a,i02b^ q Meroni^s^, 

F. S. Merritt^i, H. Merritt^^^ A. Messina^o.f, J. Metcalfe^^, A.S. Mete^^^^ C. Meyer^^, 

C. Meyer^i, J-P. Meyer^^e^ Meyer^^^, J. Meyer^^, S. MichaF, L. Micu^^^, 
R.P. Middleton^^g^ g^ Migas^^, L. Mijovic^^^, G. Mikenberg^'^^^ jyj Mikestikova^^s ^ 
M. Mikuz'^^ D.W. Miller^^, R.J. MillerS^ W.J. Millsl^^ C. Mills^^ A. Milov^'^^^ 

D. A. Milstead^^'^^'i^'^'', D. Milstein^^^^ a. A. Minaenko^^s^ Minano Moyal6^ 

LA. Minashvili'^^, A.I. Mincer^°^, B. Mindur^^, M. Mineev'^'', Y. Ming^'^^, L.M. Mir^^^ 

G. Mirabellii32a^ Mitrevski^^?^ y A. Mitsou^''^ S. Mitsui''^ P.S. Miyagawa^^Q^ 
J.U. Mj5rnmark79, T. Moa^^^a.ueb^ y_ Moeller^s, K. Monig^^^ Moser^^, 

S. Mohapatral4^ W. Mohr^^^ R. Moles- Valls^'^^ A. Molfetas^o, J. Monk'^^ E. Monnier^s, 

J. Montejo Berlingen^^^ Monticelli^° , S. Monzani^Oa.sob^ Moore^, 

G.F. Moorhead^'5, C. Mora Herrera'^'^, A. Moraes^^, N. Morange^^'^, J. Morel^^, 

G. Morello^'^'^'^^^, D. MorenoS\ M. Moreno Llacer^'^^, P. Morettini^"'", M. Morgenstern^^ , 

M. Morii^^, A.K. Morley^°, G. Mornacchi^^, J.D. Morris'''^, L. Morvaj^°\ H.G. Moser*^^, 

M. Mosidze^^^, J. Moss^°^, R. Mount^^^ E. Mountricha^°'^ S.V. Mouraviev^^'*, 

E. J.W. Moyse^'^, F. Mueller^^^^, J. Mueller^^s^ K. Mueller^^, T.A. M^ille^9^ T. Mueller^\ 
D. Muenstermann^o, Y. Munwes^^^, W.J. Murray^^g^ i Mussche^"^, E. Musto^^^^ 

A.G. Myagkov^28^ j^I. Myska^^s^ q Nackenhorst^^, J. Nadal^^^ ^ Nagai^*^°, R. Nagai^^^ 
K. Nagano'^^, A. Nagarkar^^^, Y. Nagasaka^^, M. Nagel'^^, A.M. Nairz^", Y. Nakahama^o, 
K. Nakamura'^^, T. Nakamura^^^, I. Nakano^^°, G. Nanava^^ A. Napier^'^^ 
R. Narayan^s^, M. Nash'^'^''', T. Nattermann^i , T. Naumann''^, G. Navarro^*^^^ 
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H.A. NeaF, P.Yu. Nechaeva^^, T.J. Neep^^^ A. Negri^i^^'ii^*^, G. NegrF, M. Negrini^O'^, 

S. Nektarijevic'^9, A. Nelson^ss, T.K. Nelson^'^^ S. Nemecek^^s^ Nemethyios, 

A.A. Nepomuceno^'^'', M. Nessi^^'"'', M.S. Neubauer^'^^ , M. Neumann^'^^, A. Neusiedl^^ 

R.M. Neves^°^, P. Nevski^^, F.M. Newcomer^^o^ Newman^^, V. Nguyen Thi Hong^^'^, 

R.B. Nickerson^^^, R. Nicolaidou^^^ , B. Nicquevert^° , F. Niedercorn^^^, J. Nielsen^^'', 

N. Nikiforou^s, A. Nikiforovi^, V. Nikolaenko^^s ^ Nikolic- Audit ^^ K. Nikolics''^ 

K. Nikolopoulos^^, H. Nilsen^^, R Nilsson^, Y. Ninomiya^^^ , A. Nisati^^^^, R. Nisius^'^, 

T. Nobe^^'^, L. Nodulman'^, M. Nomachi^^*^, I. Nomidis^^^, S. Norberg^^\ M. Nordberg^^, 

J. Novakova^^''', M. Nozaki^^, L. Nozka^^'^, A.-E. Nuncio-Quiroz^^ , G. Nunes Hanninger^^, 

T. Nunnemann^^ E. Nurse'^^ B.J. 0'Brien^^ D.C. O'Neil^^^^ V. 0'Shea^^ L.B. Oakes^^, 

E.G. Oakham29'/, H. Oberlack^^, J. Ocariz'^8^ A. Ochi^^^ S. Oda'^^ S. Odaka^^, 

J. Odier^^^ H. Ogren^o, A. Oh^^, S.H. Oh4^ C.C. Ohm^o, T. Ohshima^o\ 

W. Okamura^^^, H. Okawa^^, Y. Okumura^^, T. Okuyama^^^, A. Olariu^^'^, 

A.G. Olchevski^^, S.A. Olivares Pino^^a^ 01iveira^24a,i^ j3_ Qliveira Damazio^^, 

E. Oliver Garcial6^ D. Olivito^^o^ Olszewski^^, J. Olszowska^^, A. Onofre^24a,a6^ 

P.U.E. Onyisi^i''*^ C.J. Oram^^^a^ Qreglia^i, Y. Oren^^^^ D. Orestanoi34a,i34b^ 

N. Orlando^2a,72b^ q Qropeza Barrera^^^ R.S. Orr^^^, B. Osculati^o^'^o'", R. Ospanov^^^, 

C. Osuna^^ G. Otero y Garzon^^, J.R Ottersbach^o^ M. Ouchrif^^sd^ ^ Ouellette^^^ 
E. Ould-Saada^^''', A. Ouraou^^^, Q. Ouyang^^*^, A. Ovcharova^^, M. Owen^^, S. Owen^^^, 
V.E. Ozcan^^'', N. Ozturk^, A. Pacheco Pages^^, C. Padilla Aranda^^, S. Pagan Griso^^, 
E. Paganis^39^ q PahP, F. Paige^^, P. PaisS^ K. Pajchel^l^ G. Palacino^^^'^, 

CP. Paleari^, S. Palestini^o, D. Pallin^^, A. Palmai24a^ j.d. Palmeri^, Y.B. Pani'^^ 

E. Panagiotopoulou^'^ , J.G. Panduro Vazquez'^^, P. Pani^''^, N. Panikashvili^ , 
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